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Abstract 



We give a method to construct cospectral graphs for the normalized Laplacian 
by a local modification in some graphs with special structure. Namely, under some 
simple assumptions, we can replace a small bipartite graph with a cospectral mate 
without changing the spectrum of the entire graph. We also consider a related result 
for swapping out biregular bipartite graphs for the matrix A + tD. 

We produce (exponentially) large families of non-bipartite, non-regular graphs which 
are mutually cospectral, and also give an example of a graph which is cospectral with 
its complement but is not self-complementary. 

Keywords: normalized Laplacian; cospectral; bipartite subgraph swapping 



1 Introduction 

Spectral graph theory examines relationships between the structure of a graph and the 
eigenvalues (or spectrum) of a matrix associated with that graph. Different matrices are 
able to give different information, but all the common matrices have limitations. This is 
because there are graphs which have the same spectrum for a certain matrix but different 
structure — such graphs are called cospectral with respect to that matrix. 
The following are some of the matrices studied in spectral graph theory: 

*This work was done with support of an NSF Mathematical Sciences Postdoctoral Fellowship while at 
UCLA. 
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xn<c<i<^AAnA 

(a) Cospectral for A (b) Cospectral for L (c) Cospectral for Q (d) Cospectral for C 

Figure 1: Examples of cospectral graphs for A, L, Q and L. 

• The adjacency matrix A. This is defined by A(u,v) = 1 when u and v are adjacent 
and otherwise. The spectrum of the adjacency matrix can determine the number of 
edges and if a graph is bipartite, but it cannot determine if a graph is connected (see 
[6] for more information). See Figure [Taj the "Saltire pair" [5j [15], for an example of 
two graphs which are cospectral with respect to the adjacency matrix but not any of 
the other matrices we will discuss. 

• The Laplacian L = D — A. This matrix is also known as the combinatorial Laplacian 
and is found by taking the difference of the diagonal degree matrix D and the adjacency 
matrix. The spectrum of the Laplacian can determine the number of edges and the 
number of connected components, but it cannot determine if a graph is bipartite (see 
j6] for more information). See Figure [TBI [T0[ IT5] . for an example of two graphs which 
are cospectral with respect to the Laplacian but not any of the other matrices we will 
discuss. 

• The signless Laplacian Q — D + A. This matrix is found by taking the sum of the 
diagonal degree matrix D and the adjacency matrix. The spectrum of the signless 
Laplacian can determine the number of edges and the number of connected components 
which are bipartite, but it cannot determine if a graph is bipartite or connected (see 
El for more information). See Figure [Tel |15j . for an example of two graphs which 
are cospectral with respect to the signless Laplacian but not any of the other matrices 
we will discuss. 

• The normalized Laplacian C. This matrix is defined by C = D^^LD" 1 ^ 2 , where by 
convention if we have an isolated vertex then it will contribute to the spectrum. 
The spectrum of the normalized Laplacian is closely related to the spectrum of the 
probability transition matrix of a random walk. This spectrum can determine if a 
graph is bipartite and the number of connected components, but it cannot determine 
the number of edges (see [I] for more information). See Figure [Tdl for an example of 
two graphs which are cospectral with respect to the normalized Laplacian but not any 
of the other matrices we will discuss. 

One way to understand what structure the spectrum of a matrix cannot identify is to 
study cospectral graphs. Cospectral graphs for the adjacency matrix [T2J Q21 HU [T51 [TB] 
and the Laplacian matrix [T51 [T71 [20] have been studied, particularly for graphs with few 
vertices. Cospectral graphs for the signless Laplacian have been little studied beyond their 
enumeration and rules which apply to all of these matrices [T5] . 
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Little is also known about cospectral graphs for the normalized Laplacian compared to 
other matrices. Previously, the only cospectral graphs were bipartite (complete bipartite 
graphs [31 [20] and bipartite graphs found by "unfolding" a small bipartite graph in two ways 
[2]) or were regular and cospectral for A (since a regular graph is cospectral for all of these 
four matrices if it is cospectral for any one of them). 

In Table [H we have listed the number of graphs with a cospectral mate for the normalized 
Laplacian matrix for graphs on nine or fewer vertices, counted using Sage (see Appendix [A] 
for example code). The number of cospectral graphs for the normalized Laplacian have also 
been counted by Wilson and Zhu [21j, though they give percentages and not the count of 
how many have cospectral mates. We have also included similar counts for the other three 
matrices which come from [Tj [T5] . 
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Table 1: Number of graphs with a cospectral mate for the various matrices 

Given the relatively small number of graphs with a cospectral mate with respect to the 
normalized Laplacian, it is surprising that so little is known about forming cospectral graphs 
for that matrix. The problem is that some of the main tools that are used to form cospectral 
graphs for other matrices do not generalize to the normalized Laplacian. One such example 
is a technique known as switching, which is accomplished by replacing edges by non-edges 
and non-edges by edges between two subsets (see [H E51 EB] ; a simple example is shown in 
Figure [2]). Given some basic assumptions, this is an easy way to construct cospectral graphs 
for the adjacency matrix. However, switching does not in general work for the normalized 
Laplacian, in particular, it will only be guaranteed to work when the degrees are unchanged 
(see [3]). So the graphs shown in Figure [2] are not cospectral with respect to the normalized 
Laplacian. 




Figure 2: An example of two cospectral graphs for the adjacency matrix related by switching. 



3 



In Section [2] we will introduce a method to construct cospectral graphs for the normalized 
Laplacian. This construction will work similarly to switching in that we will make a small 
local change to the graph by swapping in one bipartite graph with a cospectral mate and 
show the two graphs still share the same eigenvalues. In Section Owe will show that if we add 
additional constraints to the bipartite graphs which are swapped then the resulting graphs 
are also cospectral with respect to A, L and Q. In Section H] we will show how to use this 
construction to produce large families of graphs which are mutually cospectral. Finally, in 
Section [5] we will give some concluding remarks. 

2 Swapping bipartite subgraphs 

The method of finding cospectral graphs for the adjacency matrix by switching reduces to 
making a local change of the graph from one bipartite subgraph to its complement and 
showing that the spectrum is unchanged. We will consider something similar, namely a local 
change, but instead of replacing a bipartite subgraph with its complement we will swap out 
a bipartite subgraph for a cospectral mate. 

For a subset W of the vertices V of G, we will let G[W] be the induced subgraph of G 
on the vertex set W. 

Theorem 1. Let Pi and P 2 be bipartite, cospectral graphs with respect to the normalized 
Laplacian on the vertex set B U C such that all edges go between B and C and where all 
vertices in B have degree k for both graphs. 

Let G\ be a graph on the vertices A U A' U B U C where G\ [A U A'] is an arbitrary graph; 
Gi[B] and Gi[C] have no edges; there are no edges going between A and B, between A and 
C , or between A' and C ; G\ [A' U B] is a complete bipartite graph; and G\ [B U C] = P±. The 
graph G2 is defined similarly except that G2IB U C] = Pi- Then G\ and G2 are cospectral 
with respect to the normalized Laplacian. 

If the dimension of the eigenspace associated with A = 1 intersected with the subspace of 
vectors that are nonzero only on B is the same for both Pi and P2, then the graphs Hi and 
H 2 are also cospectral with respect to the normalized Laplacian, where Hi only differs from 
Gi in that Hi [B] is the complete graph and similarly H 2 only differs from G 2 in that H 2 [B] 
is the complete graph. 

An example of the construction described in Theorem [1] is shown in Figure [3j In this 
case Pi and P 2 are the graphs K11 U K\$ and U Ki t3 where we have put the degree 
1 vertices into B. It is well known that the spectrum of a complete bipartite graph K PA 
is ()W, \\ p+q ~ 2 \ 2^, where the exponent indicates multiplicity, hence these graphs are easily 
seen to be cospectral. Moreover it is easy to check that the dimension of the eigenspace for 
the eigenvalue 1 intersected with the subspace of vectors which are nonzero only on B is 5 
for both graphs, i.e., we could also have put a complete graph on the vertices of B in the 
graphs shown in Figure [3] and still had a cospectral pair. 
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B 

(a) Gi 



B 

(b) G 2 



Figure 3: An example of cospectral graphs using the construction given in Theorem [TJ 

Additional examples of graphs Pi and P2 satisfying (both) the conditions on Theorem [1] 
are P\ = Kf . p U K^ q and P2 = Ke, r U Kg tS where the vertices of degree £ are all placed in B 
and p + q = r + s. Two further examples are shown in Figures H] and [5j 




(a) Pi (b) P 2 

Figure 4: A pair of cospectral graphs satisfying the conditions of Theorem [TJ 

To prove Theorem [TJ we will find it convenient to work with the harmonic eigenvectors of 
the normalized Laplacian. Namely, if x 7^ is an eigenvector associated with the eigenvalue 
A, i.e., £x = Ax, then the corresponding harmonic eigenvector is y = D~ x / 2 x. This translates 
the relationship £x = Ax into (D — A)y = XDy, which at a vertex v becomes 

5>(u) = (l-A)y(^), (1) 

where d(v) is the degree of the vertex v. We will say that two harmonic eigenvectors yi and 
y2 are orthogonal if y^-Dyi = 0, i.e., if the corresponding eigenvectors are orthogonal. 

Before we begin the proof of Theorem [TJ it will be useful to make some observations about 
bipartite graphs that will come up in the proof. 

Lemma 2. Let P be a bipartite graph on the vertices BUC where all edges go between B and 
C and the vertices in B all have degree k. Further, i/x is an eigenvector for the normalized 
Laplacian then we can write it as x = b + c where h is the vector x restricted to the vertices 
of B, and similarly for c. Then the following hold: 
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( a ) p i 



(b) P 2 



Figure 5: A pair of cospectral graphs satisfying the conditions of Theorem [TJ 

(a) If £x = x, i.e., x is an eigenvector for the eigenvalue 1, then Ch = b and Cc = c. In 
other words, we can divide the eigenspace for the eigenvalue 1 between those which are 
nonzero only on the vertices of B and nonzero only on the vertices of C . 

(b) If x = b + c is an eigenvector associated with \, then b — c is an eigenvector associated 
with 2 — A. 

(c) // xi = 1# + lc, . . . , Xj = bj + Cj, . . . , x m = 1 B — lc are a set of orthogonal harmonic 
eigenvectors of P, then for 1 < i < m we have bj is orthogonal to 1 B and Cj is 
orthogonal to lc- In particular, for 1 < i < m the sum of the entries of hi is 0. 

Parts (a) and (b) easily follow from ([T]). For (c) we note that Xj will be orthogonal to 
the vectors |(xx + x m ) = 1 B and |(xi — x m ) = l c , also since the vertices in B are regular 
the result on the sum of the entries of bj follows by the definition of orthogonality. 

Proof of Theorem^ Let x x = 1 B + lc> • • • , x « = b» + c», . . . , x m = 1 B — l c be a full 
set of orthogonal harmonic eigenvectors of Pi associated with Ai = 0, . . . , Aj, . . . , A m = 2 
respectively. Let s denote the degree of a vertex in B for the graph G\ and k the degree 
in P\. Then for 1 < % < m we have that 5q = bj + ^/j-Ci is a harmonic eigenvector for 



G\ associated with eigenvalue ji — 1 — (1 — Aj)y-. To see this we need to examine what 



happens for each vertex v in G\ using (CQ). 

• For v & A: Then all of the entries of the vertex and its neighbors are all and so (pQ) 
trivially holds. 

• For v G A': All of the neighbors in A and A' are while the sum of the entries in B 
will be by Lemma EJ^c), so both sides again are and ([I]) holds. 

• For We note that all of the nonzero elements of % adjacent to v are in C, and 
further we have that for the harmonic eigenvector Xj for Pi that 
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So we have for G\ that 



= (1 -li)%(v)d(v), 

so ([1]) again holds. 

For »gC: We proceed similarly as we did for vertices in B and we get 



^2%(u) = y^Xj(tt) = (1 - Ai)xj(?;)d(t)) = Y^-(l - 



)%(v)d(v) = (1 -7i)xi(t;)d(u). 



Next we note by Lemma Wis) that 1# and lc are orthogonal to each of the vectors Xj 
for 1 < i < m. On Pi, the dimension of the space orthogonal to the % for 1 < i < in 
restricted to Pi is two. A basis for the orthogonal complement of spanjxj} is {1b, lc}; 
so any harmonic vector orthogonal to all of these must be a linear combination of 1b and 
lc- In particular, if we let y be a harmonic eigenvector of G\ that is orthogonal to all of 
the Xj for 1 < i < m, then y restricted to B and C is orthogonal to the Xj, so we have 
y\buc — bis + clc for some constants b and c. In other words, the harmonic eigenvector y 
is constant on the vertices of B and constant on the vertices of C. 

Everything that we have done for G\ carries over to G*2. In particular, since P\ and P 2 
are cospectral, then the 7j found by generalizing the harmonic eigenvectors of Pi and P 2 
remain the same. Furthermore, every other harmonic eigenvector orthogonal to the % must 
be constant on the vertices of B and the vertices of C. To finish off the first claim of the 
theorem, we now only need to observe that any harmonic eigenvector which is orthogonal 
to all of the Xj in G\ is also a harmonic eigenvector for G 2 for the same eigenvalue. Again 
to see this we only need to consider what happens for each vertex. 

• For v G A U A': Then all of the entries of the vertex and its neighbors are the same for 
both graphs, so (JT|) trivially holds. 

• For ugB: All of the neighbors in A' are the same, and we have the exact same number 
of neighbors in C as before with the same value and so again ([T|) holds. 

• For v G C: For the graph G\, let (3 be the fixed value of the vertices in B. We have 

= d(tO0 = (i-A)y(«W«)- 

Of course the d(v) terms cancel and we are left with (3 = (1 — X)y(v). So even though 
the degree of the vertex in C might change, it will have no effect on this relationship 
and so ((H) again holds. 
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In summary, we were able to find m — 2 harmonic eigenvectors for G\ and G 2 that gave the 
same set of eigenvalues. For any other harmonic eigenvector orthogonal to these, the same 
harmonic eigenvector worked for both graphs and so the remaining set of eigenvalues also 
agreed. So we can conclude that G\ and G2 are cospectral with respect to the normalized 
Laplacian. 

We now turn to the second statement of the theorem. Let s be the degree of a vertex 
in B in the graph H\. So let xi = 1b + lc, • • • , x i = b< + c$, . . . , x m = 1b — lc be a full 
set of orthogonal harmonic eigenvectors of Pi associated with Ai = 0, . . . , Aj, . . . , A m = 2 
respectively. First, we consider the harmonic eigenvectors of the graph Pi in the eigenspace 
corresponding to A = 1. By Lemma |2Ja) we can assume that the vectors Xj are either of the 
form bj or c,. Let the corresponding harmonic eigenvector of H% be % = Xj, i.e., we simply 
expand Xj to be zero outside of P\. Then we have two cases. 

• If Xj = bj: In this case, we know there are no problems for the vertices in A (every 
term in ([T]) is 0) or C (since ([I]) reduces to what was done in Pi). For the vertices in 
A' we can use Lemma [21^c) to see that both sides of (PQ) are 0. Finally, suppose that v 
is a vertex in B. Then we have 

y^%{u) = y^Xi(w) -Xi(f) = -%(v) = —Xi(v)s = —%(v)d(v). 
^— ' s s 

In particular, we have that 1 — 7, = — - is an eigenvalue, i.e., that 7$ = 1 + - is an 
eigenvalue, for this harmonic eigenvector. 

• If Xj = Cf. In this case we similarly know there are no problems for vertices in A, A' 
and B. For vertices in C, reduces to what we had in Pi and so we can conclude 
that this is a harmonic eigenvector for the eigenvalue 7$ = 1. 

Now suppose that Xj = bj + Cj is associated with Xi ^ 1 for Pi. Then we now will create 
two harmonic eigenvectors for Hi, namely 



y I = hi + tiCi = bi 
= hi+ t 2 Ci = bi 



1 + y / l + 4(l-A i ) 2 sfc 

2k(l - Xi) 
1 - y/l + 4(1 - Xifsk 



2k{l-Xi) 

These will be associated with the eigenvalues of 



2s + 1- Jl + A(l- Xi) 2 sk , 2s + 1 + Jl + 4(1 - Xi) 2 sk 
7i = and 7 2 = , 

' 2s ' 2s 

respectively, for Hi. Some simple computations show that the following relationships hold: 
*i(l-7i) =t 2 (l-7 2 ) = 1-Ai, ti(l-Xi)k-l = (1-^)3, and t 2 (l - XJk - 1 = (1 - l2 )s. 
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(Before we proceed to the next step of showing that these are indeed harmonic eigenvec- 
tors and eigenvalues we first make the observation that we are not creating more harmonic 
eigenvectors than we had before. This is because in Lemma 0(b) we can pair up harmonic 
eigenvectors and eigenvalues. In particular, we would have generated the same new harmonic 
eigenvectors and eigenvalues if we had used 2 — Aj and b; — Cj. So really we are taking pairs 
of eigenvectors and eigenvalues to new pairs of eigenvectors and eigenvalues.) 

To verify that these are harmonic eigenvectors with the specified eigenvalues we need to 
examine what happens for each vertex v in H\ using ([!]). We will step through and 71, 
the arguments for yf and 72 are the same. 

• For v G A: Then all of the entries of the vertex and its neighbors are all and so ([I]) 
trivially holds. 

• For v G A': All of the neighbors in A and A' are while the sum of the entries in B 
will be by Lemma EJ^c), so both sides again are and ([1]) holds. 

• For v G B: We note that all of the nonzero elements of % adjacent to v are in B and 
C, and further we have that for the harmonic eigenvector Xj for P x that 

2jxi(w) = (1 - Xi)Xi(v)k. 

u£C 

So we have for Hi that 

^Jy-j 1 (u) = t x y^Xj(tt) + ^Jxi(u) - -Xi(v) = tx(l - Xi)xi(v)k - Xi(u) 
u~u ueB 

= (tx(l - \i)k - l)x<(i;) = (1 - 7 i)y- (v)s = (1 - 7i)yKv)d(v). 

so (DO) again holds. 

• For v G C: We proceed similarly as we did for vertices in B and we get 

$>*( M ) = E x ^ M ) = 0--^)Mv)d(v)=t 1 {l- 1 i)Mv)d(v) = (l-jx)yKv)d(v). 

The remainder of this case now proceeds as before. Namely, everything that we did for 
H\ carries over for H2 (counting multiplicity of eigenvalues). Further, any other harmonic 
eigenvector orthogonal to the ones given must be constant on B and C and so any other 
harmonic eigenvector which works for H x also works for H 2 . Therefore we can conclude that 
they have the same set of eigenvalues, i.e., are cospectral with respect to the normalized 
Laplacian. □ 
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3 Swapping biregular bipartite subgraphs 



The graphs generated by Theorem [T] will generally not give cospectral graphs for other 
matrices. So it is instructive to examine the proof and try to understand the point at which 
the fact that we were using the normalized Laplacian as compared to some other matrix 
came into play. The key is understanding how (pQ) remains true even when the degrees of 
vertices in C change — if the degree of a vertex in C changes, the d(v) term on the right side 
of ([1]) will change proportionally to the sum on the left so that the equality still holds. 

In the case when the degrees of vertices do not change, for example, when Pi and P 2 
are both (k, £)-biregular, then the proof generalizes. Recall that a bipartite graph is (k,£)- 
biregular if the vertices can be partitioned into BUC where all edges go between B and C, 
the vertices in B all have degree k, and the vertices in C have degree £. We now have the 
following theorem. 

Theorem 3. Let P\ and P 2 be bipartite, cospectral with respect to the adjacency matrix, 
(k,£)- biregular graphs on the vertex set B U C such that all vertices in B have degree k and 
all vertices in C have degree £ and edges go between B and C. Further, the dimension of 
the eigenspace associated with A = for the adjacency matrix intersected with the subspace 
of vectors which are nonzero only on B is the same for both Pi and P 2 . 

Let Gi be a graph on the vertices A U A' U B U C where Gi [A U A'} is an arbitrary graph; 
Gi[B] and Gi[C] have no edges; there are no edges going between A and B, between A and 
C, or between A' and C; Gi [A' U B] is a complete bipartite graph; and Gi [BUG] — Pi- The 
graph G 2 is defined similarly except that G 2 [B U C] — P 2 . Then Gi and G 2 are cospectral 
with respect to the matrix A + tD for t arbitrary. 

Similarly, the graphs Hi and H 2 are also cospectral with respect to A + tD for t arbitrary, 
where Hi only differs from Gi in that Hi[B] is the complete graph and similarly H 2 only 
differs from G 2 in that H 2 [B] is the complete graph. 

Note when t = then A + tD = A, when t = 1 then A + tD = A + D = Q, and when 
t = —1 then A + tD = A — D = — L, so that any such pairs of graphs are cospectral with 
respect to all of the matrices A, L, and Q. Further, these are also cospectral for C since if 
t = —A — 1, we have 

det(A + (-A - l)D) = det((A - D) - \D) 

= det(D 1/2 (C-XI)D 1/2 ) 
= (Q^)det(£-AJ). 

Since the two graphs have the same degree sequence, the product Yl di is the same, so the 
characteristic polynomials also are the same. 

In order to use this theorem, we must find two (k, £)-biregular bipartite graphs which are 
cospectral and for which the dimensions of the eigenspace of restricted to B agree on both 
graphs. There is one special case for which this is much easier, namely the theorem does 
not prohibit the possibility that Pi and P 2 are the same graph (so trivially are cospectral so 
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we only are reduced to checking the dimension of the eigenspaces) . This is not so trivial as 
it might sound since we are distinguishing the two parts of the bipartite graph in the proof 
of the theorem, so that while P\ and P2 are the same graph we are attaching them in two 
different ways and so the resulting graphs might not be isomorphic. 

As an example of this consider the graphs shown in Figure These are both the same 
graph (just flipped) and a simple check shows that the dimension of the eigenspace associated 
with and restricted to B is 2 in both graphs. Therefore Theorem [3] applies, and any two 
graphs to which these two graphs are attached are cospectral. But note that in P2 there are 
two vertices in C with the same neighbors (marked), but there are no such pairs in P\. So 
the resulting graphs are not isomorphic as long as A' 7^ 0. 




(a) A 




(b) P 2 



Figure 6: A pair of cospectral graphs satisfying the conditions of Theorem [31 



So for example if we consider the graphs Hi and H2 where A U A' = and we induce a 
complete graph on the vertices of B, then Theorem [3] shows that the resulting graphs are 
cospectral, and again since there are two vertices with degree 3 sharing common neighbors in 
one graph but not the other they are non-isomorphic. These graphs are shown in Figure [7J 
These graphs also happen to be complements of one another, and so give an example of a 
graph which is cospectral with its complement but not self-complementary. (Note the graph 
shown in Figure [Tel is another example of a graph which is cospectral with its complement, 
but in that case it is only cospectral with respect to Q and not with respect to A, L or C) 




(a) H x (b) H 2 



Figure 7: A pair of non-isomorphic cospectral graphs for A, L, Q and C (Note Hi = H 2 .) 

The proof for Theorem [3] will proceed similar to the proof for Theorem [TJ so we will skip 
some of the routine computations and provide an outline of the proof. Note that if x is an 
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eigenvalue associated with eigenvalue A, i.e., (A + tD)x = Ax, then at a vertex v we have 

5>(u) = (A-fd(«))x(t;). (2) 

We first start with some simple properties of bipartite graphs that will be helpful in the 
proof. The proofs are similar to Lemma [2] and we will omit them. 

Lemma 4. Let P be a bipartite (k,£)-biregular graph on the vertices B U C where all edges 
go between B and C. Further, z/x is an eigenvector for the adjacency matrix, then we can 
write it as x = b + c where b is the vector x restricted to the vertices of B, and similarly 
for c. Then the following hold: 

(a) If Ax = 0, i.e., x is an eigenvector for the eigenvalue 0, then Ah = and Ac = 0. In 
other words, we can divide the eigenspace for the eigenvalue between those which are 
nonzero only on the vertices of B and nonzero only on the vertices of C . 

(b) J/x = b + c is an eigenvector associated with \, then b — c is an eigenvector associated 
with —A. 

(c) 7/xi = \fkl B + vile, ■ ■ ■ , Xj = bj + Cj, . . . , x m = \/kl B — y/ile are a set of orthogonal 
eigenvectors of P, then for 1 < i < m we have bj is orthogonal to 1 B and Cj is 
orthogonal to lc- In particular, for 1 < i < m the sum of the entries of hi is 0. 

Proof of Theorem We will first consider the case when the vertices on B induce an empty 
graph. 

Let x x = Vkl B + y/Ilc, • • • , Xj = bj + Cj, . . . , x m = Vkls — V*lc be a full set of 
orthogonal eigenvectors of Pi for the adjacency matrix associated with the eigenvalues Ai = 
\/k£, . . . , Aj, . . . , —ykl respectively. Let s denote the degree of a vertex in B in G\, and note 
by construction that i will be the degree of a vertex in C in G\. We now show how to extend 
each of these eigenvectors to an eigenvector for G\. We have the following: 

• For an eigenvector of the form x = b (i.e., only nonzero on the vertices of B) and 
A = then y = b is an eigenvector for A + tD on G\ associated with the eigenvalue of 
7 = ts. 

• For an eigenvector of the form x = c (i.e., only nonzero on the vertices of C) and A = 
then y = c is an eigenvector for A + tD on G\ associated with the eigenvalue of 7 = tl. 

• For an eigenvector x = b + c associated with an eigenvalue A 7^ 0, we can construct 
two new eigenvectors, namely 

yi = 2Ab + ( - t(s -£) + v / t 2 (s-£) 2 + 4A 2 )c and 
y 2 = 2Ab + ( - t{s -£)- ^t 2 (s- £) 2 + 4A 2 )c. 
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These are associated with the eigenvalues of A + tD for G\ 



t(s + £) + Jt 2 (s-£) 2 + 4\ 2 , t(s + £) - Jt 2 (s - £) 2 + A\ 2 
7i = — and 72 = 

respectively. 

(Note, as in the proof of Theorem [TJ that we are not creating too many new eigenvectors 
and eigenvalues, i.e., by Lemma HJ^b) we can pair up the eigenvalues and eigenvectors. 
In particular, x = b + c and A will create the same new eigenvectors and eigenvalues 
as x = b — c and —A.) 

The verification that each of these is an eigenvector /eigenvalue pair for G\ reduces to 
verifying (j2J) for each vertex, similarly as was done in Theorem [TJ We will skip these routine 
computations. 

Next we note by Lemma HJc) that 1# and lc are orthogonal to each of the new eigen- 
vectors for 1 < i < m. Since the dimension of the space orthogonal to the new eigenvectors 
restricted to P\ is two, any vector orthogonal to all of these new vectors must be a linear 
combination of 1b and lc when restricted to P\. In particular, if we let y be an eigenvector 
of G\ that is orthogonal to all of the new eigenvectors, then y restricted to B and C must be 
of the form y\buc = bis + clc for some constants b and c. In other words, the eigenvector 
y is constant on the vertices of B and constant on the vertices of C. 

Everything that we have done for G\ carries over to G^- in particular, since P\ and P2 
are cospectral, then the newly found 7, are the same for both graphs, and further, every 
other eigenvector orthogonal to the ones found must be constant on the vertices of B and 
the vertices of C. To finish off the first case of the theorem, we now only need to observe that 
any eigenvector which is orthogonal to all of the new eigenvectors in G\ is also an eigenvector 
for G<i for the same eigenvalue. This again is done in the exact same way as in Theorem [TJ 
and we will skip the computations here. 

In summary, we were able to find m — 2 eigenvectors of A + tD for G\ and G2 that 
gave the same set of eigenvalues. For any other eigenvector orthogonal to these, the same 
eigenvector worked for both graphs and so the remaining set of eigenvalues also agreed. So 
we can conclude that G\ and G2 are cospectral with respect to the matrix A + tD. 

Now we turn to the case when the vertices on B induce a complete graph. 
Using the same notation as in the previous case, we will extend the eigenvectors of P\ 
for the adjacency matrix to eigenvectors of H\ for A + tD. We have the following: 

• For an eigenvector of the form x = b (i.e., only nonzero on the vertices of B) and 
A = then y = b is an eigenvector for A + tD on Hi associated with the eigenvalue of 
7 = ts — 1. 

• For an eigenvector of the form x = c (i.e., only nonzero on the vertices of C) and A = 
then y = c is an eigenvector for A + tD on H\ associated with the eigenvalue of 7 = tl. 
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• For an eigenvector x = b + c associated with an eigenvalue A 7^ 0, we can construct 
two new eigenvectors, namely 

yi = 2Ab + (1 - t(s -£) + y / (l-t(s-£)) 2 + 4A 2 )c and 
y 2 = 2Ab + (1 - t(s -£)- sj(l - t(s - £)) 2 + AX 2 )c. 
These are associated with the eigenvalues of A + tD for Hi 

t(s + £) - 1 + J \l - t(s - £)) 2 + 4A 2 
7i = g and 

t(s + £) - 1 - J (I - t(s - £)) 2 + AX 2 



respectively. 

(Note, as in the proof of Theorem HJ that we are not creating too many new eigenvectors 
and eigenvalues, i.e., by Lemma H^b) we can pair up the eigenvalues and eigenvectors. 
In particular, x = b + c and A will create the same new eigenvectors and eigenvalues 
as x = b — c and —A.) 

The verification that each of these is an eigenvector /eigenvalue pair for Hi reduces to 
verifying ([2]) for each vertex, similarly as was done in Theorem [TJ We will skip these routine 
computations. 

The remainder of this case now proceeds as before. Namely, everything that we did for 
Hi carries over for H 2 (counting multiplicity of eigenvalues). Further, any other eigenvector 
orthogonal to the ones given must be constant on B and C, so any other eigenvector for 
A + tD which works for Hi also works for H 2 . Therefore we can conclude that they have 
the same set of eigenvalues, i.e., are cospectral with respect to A + tD. □ 



4 Large families of mutually cospectral graphs 

We can use Theorem [1] to construct large families of mutually cospectral graphs for the 
normalized Laplacian that need not be bipartite nor have the same degree sequences. 

Example 1. Let n, k be positive integers such that mi + ■ • • + = n is a partition of 
n into k positive integer parts. Let FB(m 1 ,m2,...,mj,) be the graph on n + k vertices 
bi, . . . , b n , Vi, . . . , vj~, where the bi induce a complete graph, each Vi is only adjacent to exactly 
rrii vertices from bi, . . . , b n , and each bi is adjacent to exactly one of the Vi. The collection 
of all FB(mi,m2, . . . , m^) for partitions mi + m 2 + • ■ ■ + m^ = n give a family of cospectral 
graphs for the normalized Laplacian, i.e., the spectrum is completely determined for the 
normalized Laplacian by n and k. 
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FB (§1;1) FB (52;1) FB (4;3;1) FB (4;2;2) FB (3;3;2) 



Figure 8: Fuzzy ball graphs for n = 8 and k = 3. 

We have dubbed these graphs "fuzzy balls" since they consist of a complete graph (the 
bi) with some extra vertices attached (the fuzz). In Figure [BJ we have shown the cospectral 
family that corresponds to n = 8 and k = 3. 

To see that these are all cospectral, we observe by Theorem [T] that FB(mi, m 2 , . . . , rrik) is 
cospectral with FB(mi+mi— 1, m 2 , . . . , m.j_i, 1, m i+1 , . . . , m^) by using the pair of cospectral 
graphs Ki >mi U i£i, mi and iT 1)mi+mi _ 1 U K lf i where B induces a complete graph. Therefore, 
applying this idea k — 1 times, we can conclude that for each partition of n into k parts we 
have FB(mx, m^, ■ ■ ■ , m^) is cospectral with FB(n — k + 1, 1, . . . , 1). In particular, they are 
all mutually cospectral. 

Example 2. Let n, k be positive integers such that m\ + • ■ ■ + m\. = n is a partition of n 
into k positive integer parts. Let IS(mi, rri2, ■ ■ ■ , m^) be the graph on n + k + 1 vertices a, b\, 
. . ., b n , v\, . . ., Vk, where the bi and a induce a star graph with a as the central vertex, each 
Vi is only adjacent to exactly rrii vertices from b\, . . . ,b n , and each bi is adjacent to exactly 
one of the V{. The collection of all IS (mi, 777.2, • • • , n^k) for partitions ni\ + rri2 + ■ ■ ■ + mk = n 
give a family of cospectral graphs, i.e., the spectrum is completely determined by n and k. 

We have dubbed these graphs "inflated stars". In Figure [9], we have shown the cospectral 
family that corresponds to n = 8 and k = 3. 




IS ( £>■ 1; 1 ) IS(5p2;l) IS(43;1) IS(42;2) IS (33;2) 

Figure 9: Inflated star graphs for n = 8 and k = 3. 



Again, to see that these are all cospectral, we observe by Theorem [T] that the graph 
IS(mi,m2, . . . ,rrik) is cospectral with IS{mi + mi — l,m 2 , . . . , mj_i, l,m i+ i, . . . ,mfe). There- 
fore, applying this idea k — 1 times, we can conclude that for each partition of n into k parts, 
we have IS(rrii, m 2 , . . . , mj.) is cospectral with IS(n — k + 1, 1, . . . , 1). In particular, they 
are all mutually cospectral. 
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In both of these examples, we are using partitions to form large families. The number 
of partitions grows subexponentially with n, however it is possible to use Theorems [1] or [3] 
to construct large mutually cospectral families that grow exponentially with n, either with 
respect to the normalized Laplacian, or more generally with respect to A + tD. 

Theorem 5. For n large, there exists a family o/2'- n / 7 J non-isomorphic, mutually cospectral 
graphs with respect to the normalized Laplacian. 

Proof. Find a graph G on \n/7\ vertices that has a trivial automorphism group (for n > 42 
this is easy to do). For each vertex v of G, attach one of the two inflated stars shown in 
Figure \IU\ by identifying v with the vertex marked v in the inflated star. This constructs 
2\G\ = 2l n AJ graphs on 7[n/7\ < n vertices. We note that for any one of these graphs we can 




V V 



Figure 10: The two different inflated stars to add to vertices v in Theorem [51 

easily recover G and identify which attachment was done at v and so these graphs are all 
non-isomorphic. Further, by Theorem [1] changing our choice of widget at any single vertex 
does not change the spectrum and so all of the graphs are cospectral. □ 

We have not tried to optimize the construction. For example, if instead of using the 
family of inflated stars with n = 4 and k = 2 we used the family of inflated stars with n = 21 
and k = 5 (of which there are 101 such non-isomorphic members), then we get a family with 
rate of growth 101 L«./26J w ii^T which is faster than 2 L™/7j w 1-10 409 n . 

We similarly have the following result for graphs cospectral with respect to A + tD. 

Theorem 6. Forn large, there exists a family o/ 2 L™/ 13 J non-isomorphic, mutually cospectral 
graphs with respect to the matrix A + tD. 

Proof. Find a graph G on |_rz. / 13J vertices that has a trivial automorphism group (for n > 78 
this is easy to do). For each vertex v of G attach one of the two graphs shown in Figure [TTT 
which are taken from Figure El This constructs 2' G ' = 2 L"/ 13 J graphs on 13[n/13j < n 
vertices. We note that for any one of these graphs we can easily recover G and identify which 
attachment was done at v, so these graphs are all non-isomorphic. Further, by Theorem [3] 
changing our choice of widget at any single vertex does not change the spectrum and so all 
of the graphs are cospectral. □ 
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V 



V 



Figure 11: The two different widgets to add to vertices v in Theorem |6j 

5 Concluding remarks 

The method of swapping bipartite subgraphs was discovered by examining all cospectral 
pairs of graphs on at most eight vertices (all cospectral pairs on at most seven vertices has 
previously appeared in Tan [20]), which included several simple examples of this type. The 
set of all cospectral graphs on at most eight vertices can be easily constructed in sage (see 
Appendix |A|) . 

By looking at the small cases, we also discovered some new graphs which are regular, 
but are cospectral for the normalized Laplacian with graphs which are not regular. We have 
already seen that the four-cycle and the eight-cycle are members of cospectral pairs (see 
Figures ITdl and Also, K n>n is cospectral with K Ptq for p + q = In. Another example of 
a graph which is regular, but is cospectral with a graph which is not regular, is shown in 
Figure [121 These examples show that the normalized Laplacian cannot, in general, detect 
whether a graph is regular. 




Figure 12: A regular non-bipartite graph which is cospectral for the normalized Laplacian 
with a graph which is not regular. 

As noted in the introduction, the number of edges in a graph is not uniquely determined 
by the spectrum of the normalized Laplacian. As a consequence, it is possible for a graph 
to be cospectral with one of its subgraphs; examples of this are shown in Figure [13j Further 
examples of graphs which are cospectral with respect to the normalized Laplacian but have 
differing number of edges are shown in Figure [TH (there is currently no known method 
to generate cospectral graphs with differing number of edges other than using complete 
bipartite graphs). This also indicates the difficulty in counting the number of cospectral 
graphs. While for the adjacency and the combinatorial Laplacian, we could first subdivide 
the graphs according to the number of edges and work with this coarsening when finding 
cospectral pairs (see [15]), this is no longer possible for the normalized Laplacian. 

It would be interesting to find techniques for constructing cospectral graphs for the nor- 
malized Laplacian which differ in the number of edges. Similarly, it would be interesting to 
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(a) Cospectral pair 



(b) Cospectral pair 





(c) Cospectral pair 



(d) Cospectral pair 



Figure 13: Examples of graphs which are cospectral with a subgraph for the normalized 
Laplacian. 



investigate the cospectral graphs for the signless Laplacian and see if there are any construc- 
tions unique to forming cospectral graphs for that matrix (similar to what Theorem [1] is for 
the normalized Laplacian). 
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A SAGE code to generate cospectral graphs 

The code listing below shows how to use Sage version 4.5.2 [19] to generate and check 
cospectral graphs with 8 vertices. In the listing below, we specifically generate cospectral 
graphs with respect to the normalized Laplacian. The lines starting with the number sign 
are comment lines. 

def DinverseA (g) : 

"Calculate D~{-1}A for a graph g" 
A=g.adjacency_matrix() . change_ring(QQ) 
for i in range (g . order ()) : 

A. rescale_row(i , 1/ len ( A . nonzero_positions_in_row(i)) ) 
return A 

# Calculate the graphs and store them in cospectral_list . 

# Vie use the DinverseA function defined above to determine 

# if two graphs are cospectral. 

# This command takes a few minutes to complete . 
cospectral_list=graphs . cospectral_graphs(8, 

graphs=lambda g: min (g . degree ()) >0) 

# Give a list of all graph6 strings of these graphs. 
graph6_list= [ [g . graph6_str ing () for g in glist] 

for glist in cospectral_list] 

# Show all of the graphs, each row being cospectral. 
for glist in cospectral_list : show(glist) 
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# Get the first two cospectral graphs and check 

# to make sure they are cospectral. 
graphl , graph2 = cospectral_list [0] 

graphl_poly=graphl . laplacian_matrix(normalized=True) . charpolyO 
graph2_poly=graph2 . laplacian_matrix(normalized=True) . charpolyO 
graphl_poly==graph2_poly 

# Construct fuzzy ball graphs corresponding to the partitions 

# 4=3+1 and 4=2+2 and check that they are cospectral. 
f bl = graphs . Fuzzy B all Graph ([3,1] ,0) 

fb2=graphs . FuzzyBallGraph ( [2 ,2] ,0) 

fbl_poly = fbl . laplacian_matrix(normalized=True) . charpolyO 
fb2_poly=fb2 . laplacian_matrix(normalized=True) . char poly () 
fbl_poly==fb2_poly 
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